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We show how to use two-mode squeezed light to exponentially enhance cavity-based dispersive 
qubit measurement. Our scheme enables true Heisenberg-limited scaling of the measurement, and 
crucially, is not restricted to small dispersive couplings or unrealistically long measurement times. It 
involves coupling a qubit dispersively to two cavities, and making use of a symmetry in the dynamics 
of joint cavity quadratures (a so-called quantum-mechanics-free subsystem). We discuss the basic 
scaling of the scheme and its robustness against imperfections, as well as a realistic implementation 
in circuit quantum electrodynamics. 
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Introduction. - Research in quantum metrology has es¬ 
tablished that squeezed light and entanglement are key 
resources needed to approach truly fundamental quan¬ 
tum bounds on measurement sensitivity [T] . Perhaps the 
best known application is interferometry: by injecting 
squeezed light into the dark port of an interferometer, one 
dramatically enhances its sensitivity to small phase shifts 
[a [3], reducing the imprecision below the shot-noise limit. 
Many of these ideas for squeezing-enhanced measure¬ 
ment were first motivated by gravitational wave detection 
Hi], and have recently been implemented in current- 
generation detectors Hi- More generally, squeezed light 
has been used to enhance the measurement sensitivity in 
optomechanics [9] and even biology [10]. 

Ultrasensitive detection is also essential for quantum 
information processing where fast, high-fidelity qubit 
readout is required to achieve fault-tolerant quantum 
computation HU. A ubiquitous yet powerful approach 
is dispersive readout, where a qubit couples to a cav¬ 
ity such that the cavity frequency depends on the qubit 
state, see e.g. [12] • The readout consists in driving the 
initially empty cavity with a coherent tone, resulting in 
a qubit-state dependent cavity field which is displaced 
in phase space from the origin [see Fig. Sa)]. High- 
fidelity readout can then be obtained by measuring the 
output field quadratures. This is the standard approach 
used in state-of-the-art experiments with superconduct¬ 
ing qubits, e.g. HMl]. 

As with interferometry, one might expect that dis¬ 
persive qubit measurement could be enhanced by using 
squeezed light. The most obvious approach would be to 
squeeze the phase quadrature of the incident light [i.e. Y 
in Fig.j^b)], thus reducing the overlap between the two 
pointer states. As discussed recently in Ref. [THl the situ¬ 
ation is not so simple, as the dispersive interaction will 
lead to a qubit-dependent rotation of the squeezing axis. 
Unlike standard interferometry, this rotation is a prob- 
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FIG. 1. Phase-space representation of dispersive qubit read¬ 
out for different input states: (a) coherent state, (b) single¬ 
mode phase-squeezed state, (c) amplitude-squeezed state, 
(d) two-mode squeezed state in the QMFS (X_,Y+). The 
purple dashed lines represent the input state and the blobs 
represent the output fields. The input state is displaced along 
the X-axis and the signal is encoded in the quadrature corre¬ 
sponding to the Y -axis with homodyne detection; as depicted 
in the leftmost panel, the readout error corresponds to the 
overlap of the two marginals. Dispersive interaction with the 
qubit rotates the output field by the angle (pqb for the ground 
state |0) (in blue) and —yqb for the excited state |1) (in red). 
Ideally, one wants the output state to be phase squeezed re¬ 
gardless of qubit state [dotted “desiderata” states in (b)]; this 
is not possible when using single-mode squeezing due to the 
qubit-induced rotation. Our new QMFS scheme [panel (d)] 
does not suffer from this problem. 


lem, as optimal dispersive qubit readout involves large 
couplings and hence large rotations. Further complexity 
arises from the fact that this rotation is frequency depen¬ 
dent. The upshot is that measurement always sees the 
amplified noise associated with the antisqueezed quadra¬ 
ture of the incident light, limiting the fidelity improve¬ 
ment from using squeezing to modest values and prevent¬ 
ing true Heisenberg scaling |16j . 

Despite the above difficulties, we show in this Letter 
that it is indeed possible to substantially improve dis- 
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FIG. 2. Circuit QED implementation of QMFS dispersive 
qubit readout with two-mode squeezed states produced by a 
nondegenerate paramp (NDPA). The signal is encoded in the 
joint quadrature y+,out oc yi,out + iz.out; see text for details. 

persive qubit measurements using squeezed input states. 
Our proposed scheme involves using two-mode squeezed 
states in a two-cavity-plus-qubit system (see Fig. [^, 
which can lead to exponential enhancement of the signal- 
to-noise ratio (SNR) in dispersive measurement and 
achieves true Heisenberg-limited scaling. This is possi¬ 
ble even for large qubit-induced phase shifts, and is thus 
in stark contrast to previous schemes using two-mode 
squeezing for interferometry or qubit readout US- 
The key to our scheme is the use of a special dynami¬ 
cal symmetry, whereby two commuting collective quadra¬ 
tures exhibit a simple rotation as a function of time. As 
these quadratures commute, they constitute a so-called 
“quantum-mechanics-free subsystem” (QMFS) |T7j and 
can both be simultaneously squeezed. The upshot is 
that one can effectively make a dispersive qubit measure¬ 
ment where now the uncertainties associated with the two 
pointer states are not limited by the uncertainty principle 
[see Fig.j^d)]. Though the scheme is extremely general, 
for concreteness we explicitly discuss an implementation 
in circuit quantum electrodynamics (QED) using a trans- 
mon qubit [m, as depicted in Fig. [| 

The dynamical symmetry used in our two-mode 
scheme crucially relies on one of the cavities having an 
effective negative frequency; it is thus related to an idea 
first discussed in the context of measurement by Tsang 
and Caves m and Wasilewski et al. [^, and which has 
since been applied to other systems [minilll]. While 
many applications use the idea to suppress the effects of 
backaction [20H22] , we instead use it as an effective means 
to exploit squeezed input light. Unlike previous studies, 
we calculate here the scaling of the resulting measure¬ 
ment sensitivity, showing that one obtains Heisenberg- 
limited scaling with incident photon number. 

Dispersive measurement and standard squeezing.- We 
start by reviewing the simplest setup where a qubit dis- 
persively couples to a single-sided cavity (frequency wi) 
with the Hamiltonian H = (wi-I-xdz) a^d |12j . Stan¬ 
dard dispersive readout involves driving the input port 
of the cavity with a coherent tone at the cavity frequency 
(photon flux hoK/4, k is the cavity damping rate). As 
illustrated in Fig. [^a), as a consequence of the disper¬ 


sive coupling, the output field is rotated by the angle 
(/ 3 qb = 2 arctan( 2 x/K) if the qubit is in the ground state 
jo) and by for the excited state |1). Writing the 

output field as dout(i) = e“*“’i‘(Xout + *Fout)/ 2 , for a 
displacement along the real axis Aiout, the signal of the 
qubit state is encoded in the phase quadrature Fl>ut; this 
quadrature is then recorded with homodyne detection. 

Measuring for an integration time r corresponds 
to evaluating the dimensionless measurement operator 
M = y/k dflljut (t) ■ The signal is the qubit-state 
dependent expectation value Ms = (M) and is the 
same for all the injected states depicted in Fig. 
The imprecision noise is the variance of the noise op¬ 
erator Mat = M — Ms- The signal-to-noise ratio 
SNR = iMgjo) - M 5 J 1 ) I/((M^^ip^) -I- (M^j^^))^/^ is, 
for this coherent state dispersive readout, SNRQ(r) ~ 

I sin(/jqb|\/2hoAtr [531 [21]. As expected, the SNR is max¬ 
imized for a phase (pqh = 7 ’‘/ 2 ; it also scales as ^/nQ, akin 
to standard quantum-limit scaling in interferometry [ 1 |. 

Next, consider what happens if we instead inject a dis¬ 
placed squeezed state (squeeze parameter r) into the cav¬ 
ity. As already discussed, this is not as beneficial as one 
would hope, as one always sees the noise of the anti¬ 
squeezed quadrature (oc e^”) [T6l[25]. Consider the op¬ 
timal case tpqb = 7 r /2 which maximizes the signal. For 
large r, the noise behaves as 

(M^) ~ Kr[sin^(0)e“^’' -I- cos^(0)e^’'] 

+ 2-\/2sinh(2r) cos(20 — 37 r/ 4 ), (1) 

where we have dropped terms that decay exponentially 
with KT. The first line of Eq. 0 dominates in the long¬ 
time limit, and represents the contribution from zero- 
frequency noise in the output field. For this line, the 
choice 0 = tt/2 cancels the contribution from the ampli¬ 
fied quadrature, and leads to an exponential reduction 
in the noise compared to a coherent state drive |16| . In 
contrast, the second line of Eq. 0 describes the con¬ 
tribution from initial short-time fluctuations; the noise 
from the antisqueezed quadrature here remains even if 
0 = 7 r/ 2 . As a result, increasing r indefinitely does not 
improve the SNR; for a given t there is an optimal value 
[see Fig. (b-c)]. This then leads to generally mod¬ 
est enhancement of SNR compared to a simple coherent 
state drive m-, in particular, there is almost no improve¬ 
ment in the most relevant case where r ~ 1 /k [shaded 
region in Fig. [^a)]. Optimized squeezing leads at best 
to the scaling with input photon number, similar 

to a Mach-Zehnder interferometer driven with squeezed 

light [11 [51155]. 

Negative freqencies and two-mode squeezing. - To avoid 
having the measurement corrupted by the antisqueezed 
quadrature, one ideally wants to squeeze both quadra¬ 
tures of the input light. While this is impossible with a 
single cavity, it becomes conceivable using joint quadra¬ 
tures of two cavities. If % = {Xj iYj)l2 (j = 1,2) 
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FIG. 3. (a) SNR as a function of integration time r for differ¬ 
ent protocols: coherent state drive (black), displaced single¬ 
mode squeeze state (blue), two-mode squeezed QMFS setup 
(red). We assume an optimal dispersive shift x = ^/2; for 
the QMFS setup the cavity is presqueezed (to <C —1/k) with 
squeeze strength = 100. The coherent drive is turned on 
at r = 0. For the single-mode case, for each r we optimize 
the squeeze strength e^’’ € [1,100] and angle (see [25]). The 
QMFS scheme gives an exponential SNR enhancement, espe¬ 
cially in the most interesting regime where r ~ 1 /k (shaded 
region), (b) Integration time r required to achieve a fidelity 
F = 99.99 %, as a function of e^’’; parameters as in (a), except 
that no = 100. Black lines correspond to an unsqueezed drive, 
where the drive strength is increased such that the intracav¬ 
ity photon number is the same as in the QMFS scheme, i.e. 
no ^ no + 4sinh^ r. The solid curves correspond to the case 
of no photon losses (efficiency 77 = 1 ), while the dashed curves 
correspond to 77 = 0.9. (c) Total intracavity photon number 
needed to achieve F = 99.99 % in a measurement time r. 
Even with nonzero photon losses, the use of squeezing can 
dramatically reduce the number of intracavity photons. 

are the annihilation operators for the two cavities (in 
an interaction picture with respect to the free cavity 
Hamiltonians), we define X± = (Xi ± X- 2 )l'/^i = 

(Yi ± Y'2)l\f2. Since X- and F+ commute, they can 
be squeezed simultaneously, resulting in a two-mode 
squeezed state |26j . The relevant nonzero input-field 
noise correlators are {Xz^{t)Xz^(t')) = {Y±(t)Y±{t')) = 
F^'^''5(t — t'). We stress that such states have already 
been produced in circuit QED EllllH]. 

This squeezing by itself is not enough: we also need 
the dynamics of these joint quadratures to mimic the 
behavior of X and T in a single cavity, such that the 
two qubit states still give rise to a simple rotation of the 
vector formed by (X_ ,Y+). Such a dynamics is generated 
by the simple Hamiltonian [191 ED] 

H = \x{X+X_ + Y+Y.)a, = xi&lai - 4a2)d,. (2) 

The qubit thus needs to couple dispersively to both 
cavities, with equal-magnitude but opposite-signed cou¬ 
plings. The resulting dynamics is illustrated in Fig. 0 d): 
an incident field with ( 1 +) = 0 , (X_) 7 ^ 0 is rotated in 
a qubit-state dependent manner, resulting in an output 


field with (4-) 7 ^ 0 (i.e. the measurement signal). Note 
that the squeezed quadratures , Y^ are never mixed 
with the antisqueezed quadratures , Y- , hence this 
amplification will not limit our scheme. We also stress 
that the two cavities need not have the same frequency. 

The measurement protocol involves first turning on the 
vacuum two-mode squeezed drive at a time t = to < 
0 , and then turning on the coherent cavity drive(s) at 
t = 0. This coherent drive (which displaces along X_ 
but not Y^) could be realized by driving one or both the 
cavities. We take the optimal case where both cavities are 
driven and let ho At /8 denote the photon flux incident on 
each cavity due to the coherent drives. The measurement 
signal in y/_ can be constructed from the quadratures 
kj.out of the output field leaving each cavity. In what 
follows, we consider the limit Ktg ^ — 1, such that the 
measurement is not corrupted by any initial nonsqueezed 
vacuum in the cavity [25]. 

The measurement operator is now M = 
Jq dtY^^out{t)- As expected, one finds that this 
output quadrature is always squeezed, and hence the 
imprecision noise is always described by (M^) = 
independent of x- As desired, the noise is now expo¬ 
nentially reduced with respect to standard dispersive 
readout, leading to an exponential improvement of 
SNR, i.e. SNRr(T) = e’'SNRQ(T) for all integration 
times r. This is in stark contrast to the single-mode 
approach, where such an enhancement was only possible 
at extremely long times, kt > [cf. Eq. Q]. The SNR 
is plotted in Fig.|^a) as a function of integration time r, 
with comparisons against the single-mode squeezing 
and no-squeezing cases; our two-mode scheme realizes 
dramatic improvements in the most interesting regime 
where r is not much larger than 1/k. The integration 
time r required to achieve a measurement fidelity 
F = 1 — erfc(SNR/2)/2 of 99.99% is plotted against 
squeezing strength in Fig. db). Again, the QMFS 
scheme results in dramatic improvements. 

Heisenberg-limited scaling.- We now show that the 
SNR scales as the number of photons N used for the mea¬ 
surement rather than its square root y/N as it is the case 
for the standard dispersive readout [T]. For this, we de¬ 
fine the temporal mode H ^ /J" dt[din,l(t)+din, 2 (t)] [3] 

where the operator dinj describes fluctuations in the 
resonator-j input field. The total number of input pho¬ 
tons N = Ns Nd has a contribution from squeezing 
Ns = (A^A) = 2 sinh^r and Nd from the coherent dis¬ 
placement. Focusing on times r I/k, we can ignore 
the transient response to the coherent drive, and hence 
Nd = |hoKr. Fixing N and taking to <C —1/k, the opti¬ 
mal SNR is obtained for Ns = N'^/[2(N -j- 1)], and is 

SNRopt = 21 sin ip^t, \ N a/1 -|- 2/N 21 sin (/?qb | A^, (3) 

where we have taken the large N limit. Eq. <§ corre¬ 
sponds to true Heisenberg scaling for any value of the 
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dispersive coupling. Such scaling is not possible using 
single-mode squeezed input light (see [25]). 

Our QMFS scheme also shows an improved, 
Heisenberg-like scaling of the SNR with the intracavity 
photon number h. Note that the SNR for the QMFS 
scheme has the same form as the SNR for a standard 
(r = 0) dispersive readout made using a larger drive 
flux n^e^'". If we fix the intracavity photon number 
h = ho cos^(i^qb/2) -I- 2sinh^r and optimize r, the re¬ 
sulting SNR scales as SNRopt — 2| sin((/jqb/2)|h-\/KT, as 
opposed to the conventional SNRq, oc v^- 

Robustness against imperfections.- Our discussion of 
the QMFS scheme so far has assumed a broadband, 
pure squeezing source. The purity of the squeezing is, 
however, not crucial; our scheme is insensitive to the 
antisqueezed quadratures, and hence it is not essential 
that their variances be as small as possible. For a fi¬ 
nite squeezing bandwidth F, the input squeezing spec¬ 
trum will typically have a Lorentzian line shape |29j . We 
find that the effects of a finite bandwidth are equiva¬ 
lent to an effective reduction of the squeezing strength; 
the SNR for the scheme is simply reduced by a prefac¬ 
tor ^JTT/\^T -\- (e^’' — 1)(1 — 6“^"^)] [2S|. One thus only 
needs a modest bandwidth, e.g. F ^ 10 k is enough for 
KT ^ 10 and ~ 10. 

The lack of any enhanced Purcell decay is also cru¬ 
cial, as in our protocol the squeezing is turned on well 
before the coherent measurement tone. Having a finite 
squeezing bandwidth can in fact be an advantage as it 
helps suppress Purcell decay of the qubit. This decay 
corresponds to relaxation of the qubit by photon emis¬ 
sion from the cavity [30] ■ As typical detunings A k, 
there is a wide range of ideal squeezing bandwidths sat¬ 
isfying K <C F ^ A. Such bandwidths are large enough 
to allow a full enhancement of the SNR (with r > 1 /k), 
and small enough that the squeezing does not apprecia¬ 
bly modify cavity-induced Purcell decay (see 125]). 

Another nonideality is asymmetry in the system pa¬ 
rameters. While the two cavity frequencies can differ, 
we have assumed so far that they have identical damp¬ 
ing rates (ki = K 2 = k) and that the dispersive cou¬ 
pling strengths satisfy xi = ~X 2 = X- Deviation from 
either of these conditions breaks the symmetry yield¬ 
ing a QMFS, causing an unwanted coupling between 
the squeezed quadratures {X-,Y^) and the antisqueezed 
quadratures (A+,F_). The structure of the QMFS can 
persist in the presence of asymmetries for long measure¬ 
ment times KT 3> I, under the condition |25j 

Xl +X2 ^ - K2 

Xl - X2 Ki -H K2 ' 

The SNR enhancement can however be preserved for 
measurement times r ^ 1 /k by optimizing Sk/6x, as 
illustrated in Fig. Qa). Although this might not be 
necessary in practice, all parameters in Eq. Q can be 



FIG. 4. (a) SNR enhancement as a function of the disper¬ 

sive shift asymmetry (xi ,2 = 5x A x) for different resonator 
linewidth asymmetries (ki ,2 = k±5k) calculated for x = k/2, 
KT = 10 and e^’' = 100. The dashed line is the maximal 
SNR obtained by optimizing 5k. (b) Calculated dispersive 
shifts as a function of transmon anharmonicity Ec from a 
numerical diagonalization of transmon-resonator system for 
each of the resonators. The parameters are Ej/h = 25 GHz, 
— 7.6 GHz, ai2l2-K = 7.9GHz, — 8MHz and 

= 15 MHz. The vertical dashed line shows a typical 
value of Ec that leads to equal and opposite dispersive shifts. 


tuned in situ nainiisi] thereby greatly relaxing the 
constraints on the system. 

Finally, like any scheme employing squeezing, photon 
losses effectively replace squeezed fluctuations with ordi¬ 
nary vacuum, causing the SNR improvement to saturate 
as a function of squeezing strength [25]. Despite this, our 
scheme still yields considerable advantages for finite loss 
rates, see Figs.j^b) and (c). 

Implementation in circuit QED.- We now turn to a 
possible realization of this protocol in circuit QED. All 
parameters discussed here are readily achievable exper¬ 
imentally. As illustrated in Fig. [^ a transmon qubit 
is coupled to two resonators, one in the usual disper¬ 
sive regime (A > Ec) while the other in the “strad¬ 
dling” regime (A < Ec) [UlISS]. Here, A is the qubit- 
resonators detuning and Ec the transmon anharmonic¬ 
ity. This yields dispersive couplings x having opposite 
signs as required, see Fig. [^b). An alternative strategy 
is to use a fluxonium or a flux qubit which exhibits a 
richer dispersive shift profile [31]. Note that either ap¬ 
proach does not entail a sacrifice of qubit coherence via 
enhanced Purcell decay [25]. The displaced two-mode 
squeezed state required at the input can either be gener¬ 
ated by a NDPA such as the Josephson parametric con¬ 
verter m, a Josephson paramp m or the Bose-Hubbard 
dimer [28] . 

Conclusion.- We have presented a realistic measure¬ 
ment protocol that allows one to exponentially enhance 
dispersive measurement using two-mode squeezed light, 
enabling Heisenberg-limited scaling even with large dis¬ 
persive couplings. Our scheme crucially makes use of 
a special symmetry in the dynamics of joint cavity 
quadratures, a so-called “quantum-mechanics-free sub¬ 
system”. It could be straightforwardly generalized to 
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allow Heisenberg-limited scaling in any interferometric 
setup having large signal phase shifts. 
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I. QUBIT READOUT WITH QMFS AND TMSS 

In this section we detail the derivation of the results presented in the paper for qubit readout in a quantum- 
mechanics-free subsystem (QMFS) and driven by two-mode squeezed states (TMSS). 


A. Noise 


We consider a two-cavity setup where a single qubit is dispersively coupled to both cavities with opposite dispersive 
shifts, i.e. xi = ~X 2 = X- The cavities can have different frequencies. We suppose they have the same damping rate, 
i.e. Ki = K 2 = K- The cavities are coupled to zero-temperature environments at times t < to and, at times t > to, 
a TMSS is injected at resonance with the cavities. In the rotating frame of the cavity frequencies, we define four 
quadratures for the TMSS, X±,Y±, characterized by the correlations 
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where Z = (A_ F+ X+ F_) 

TMSS is composed of two coupled QMSS: the squeezed subsystem Z 5 = (A_ F+)^ and the anti-squeezed subsystem 
7ia = (A_|_ Y-)^. In order to involve only the squeezed QMFS in the dynamics, the Hamiltonian of the system is 
chosen to be H = |x(A_|_A_ -|- Y+Y-)<7z- As explained in the paper, this is equivalent to having opposite dispersive 
shifts. The Langevin equations for the intracavity quadratures are found from input-output theory [SI] 
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The squeezed and anti-squeezed subsystems are indeed clearly decoupled. The output quadratures are calculated from 
integrating Eq. (S2) and using the input-output boundary condition Zout = -\/kZ -|- Zin [Slj. The two-time output 
correlations of the squeezed QMFS are then equal to 
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where 1 is the 2 x 2 identity matrix and R[(/)] = I rotation matrix. Here t,t' > to and to < 0 . The 
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correlations of the anti-squeezed QMFS are found by replacing e“^’' by in Eq. (S3). The noise associated with the 
measurement operator M(r) = /J" dtl))ut,-i -(0 then reads 
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The measurement noise is identical for the two qubit states, {a^) = ±1, for all times t and to- The last terms in 
Eqs. (S3)-(S4) are due to the leakage of the initial vacuum fluctuations. As expected, the anti-squeezed quadratures 
do not contribute to the measurement noise. When the cavities are sufficiently pre-squeezed, that is when the TMSS 
is injected well before the measurement (to — 1 /k), the noise is exponentially reduced. 


B. Asymmetries 


To understand the impact of realistic asymmetries in the dispersive shifts, i.e. X 2 ~Xi; ^nd damping rates, i.e. 
7 ^ Hi, we define x = (Xi ~ X2)/2, h = {ki K 2 )/ 2 , 5x = (xi + X2)/2, 5k = (ki — K2)I‘2- In the quadrature basis, 
the Hamiltonian H = {xia\ai -f X25\a2)(7z reads 

H = ^xiX+X_ + Y+Y_)a, + lSx(Xl + (S5) 

Including the asymmetries, the Langevin equations are 
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The resulting expression for the noise in the measurement operator M{t) = y/ k fj dtYout,+ (t) is quite cumbersome. 
One can check that the measurement noise does not depend on the qubit state. The leading term for r ^ 1/ k is 
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To eliminate the contribution from the anti-squeezed QMFS, the asymmetries have to be tuned such that ^ 
Under this condition, the fnll expression for the noise is 
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The remaining terms in are constant or decaying exponentially. For finite integration times t ~ 1 /k, the 

SNR can be maximized numerically by optimizing Sx and 6k. 


C. Squeezing bandwidth 


For a finite squeezing bandwidth F, the input held correlations are 
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The noise then reads 
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The infinite bandwidth result Eq. (S4) is recovered in the limit Tr —)■ oo. 


II. QUBIT READOUT WITH SINGLE-MODE SQUEEZED STATES (SMSS) 

A. Noise 


We consider a cavity dispersively coupled to a qubit. As described for the QMFS protocol, the cavity is coupled to 
a zero-temperature environment at times t < tg and, at times t > Iq, a, SMSS is injected at resonance with the cavity. 
We note 6 the squeeze angle. In the rotating frame of the cavity frequency, we define two quadratures for the SMSS, 
A, Y, characterized by the correlations 
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where Z = {X Y)"^. The Langevin equations for the intracavity quadratures are 
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In the Fourier space and in the frame of the homodyne angle <^hj the output quadratures Zout.i^h = R'l'/^hJZout are 
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The rotation angle is frequency-dependent, the squeezed and anti-squeezed quadratures are thus simultaneously 
involved in the quadratures of the output field for any value of (/?h- For the qubit state |I), the noise of the measurement 
operator Mm{t) = ^/k /J" dthjjut.c^Ph (t) is then equal to 

(M^ ^^(r)) = KT {cos^(6* -iph- V3qb)e“^’' -bsin^(6» - iph - (/?qb)e^''} 

-b 4 sinh 2r sin (^qb cos |sin[2(6» - v^h) - fv^qb] - sin[2(6» - (/?h) - ftb’qb - X'^]} 
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— 8sinh2rcos^ ^(b>qb sin[2(0 — i^h) — fv®qb — xF ~ 2io)] sin xt sinh (S14) 

The noise {M]q |o)('^)) qubit state |0) is obtained by replacing (^qb —>■ —t/^qb and x —> —X ia Eq. (SI4). The 

anti-squeezed quadrature cannot be eliminated in the measurement noise. As a consequence, for a given measurement 
time T, there is an optimal value of the squeeze parameter r that minimizes the measurement noise. 


B. Optimal SNR 


To calculate the SNR, defined by SNR = \Ms^\o) — Mgji) |/((M^ |p^) -b the signal can be readily 

calculated from the Langevin equation Eq. (S12), giving 


^S.IO) - ^S.|i> = 2v^cos(v3h)sin(<,£)qb)«;T[I - F(r)], 


Fir) 


nr K 


sm((^qb) 


sin((^qb-b XT)e 


(S15) 



4 


The noise contribution is calculated from Eq. (S14), we take the limit to ~oo for simplicity, 
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The corresponding SNR is optimized with i^b = 0; 0 = 0 if G(t) > 0 and 6 * = | if G(t) < 0; ropt = 5 arctanh|G(T)|, 


SNRsMSs(?’opt) = V^noKT\ sin((pqb)| = V^uoktI sin(y)qb)|e’^°p* |G|r)^]i/2 ' 

The optimal squeeze strength is rather small, leading to a modest improvement of the SNR. For instance, for 
X = at T = I/k, ~ 1.8 and SNRsMSs(^opt)/SNRsMSs(0) — l-lj at r = 10 /k, ~ 3.0 and 

SNRsMSs(ropt)/SNRsMSs(0) cs 1.3; at r = 100 /k, cs 9.9 and SNRsMSs(ropt)/SNRsMSs(0) cs 2.2. 


C. Scaling 


We now focus on the situation that maximizes the signal, i.e. x = k/2. For kIq —)■ —00 and neglecting the transient 
dynamics due to the drive as well as the terms exponentially decaying with kt, the SNR reads 
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where we used = sinh^ r ~ for large squeeze parameters and N = Nd + Ng with Nd = n^KT j3,. Unlike the 
QMFS scheme, the SNR cannot be written only in terms of N and iV^; the parameter kt or hg is also involved. 

For a fixed drive strength ho, we calculate the optimal value of Ng. The SNR then scales as 

SNRsMss.opt - 81V3/Vhy^ Ng^opt = (S20) 


where we have taken the large N limit. Such a scaling, SNR oc is obtained in a Mach-Zehnder interferometer 

driven with squeezed light [S2, S3]. 

For a fixed integration time kt, we calculate the optimal value of Ng. The SNR then scales as 


SNRsMSS.opt = 1/ \/ (16A^2 _|_ i^r)KT ~ 2 Vn{kt)'^^‘’^, iV^^opt = ^ (S 21 ) 

where we have taken the large N limit. Such a scaling, SNR oc is obtained with a coherent drive without 

squeezed light [S3]. 


III. PURCELL DECAY 
A. Purcell decay due to cavity squeezing 

The coupling Hamiltonian between a cavity mode and the environment is i/diss = dw k{uj) [5'!' (w) + )>(w)] [a’^ + a], 
where a is the annihilation operator of the cavity and b{uj) the annihilation operator of the bath at oj. We apply the 
dispersive transformation to diagonalize the Jaynes-Cummings interaction between the cavity and the qubit in the 
dispersive regime, where the detuning A = Wa — Wi between the qubit and the cavity is much larger than their coupling 
strength g. Under this unitary transformation (noted with the superscript D) and at lowest order in g/A, the cavity 
operator acquires a contribution from the qubit: = a + ^d_. In the dressed basis, the qubit becomes coupled 

to the environment of the cavity. Following Ref. [S4], we apply the rotating-wave approximation in the interaction 
picture to obtain the coupling Hamiltonian between the system and the environment in the dressed basis 

Jo ^ Jo 
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The last term of the coupling Hamiltonian Eq. (S22) results in a decay channel for the qubit. The corresponding 
damping rate is 7 ^ = K{uja)- For a squeezed environnement of squeeze parameter r at resonance with the cavity 
frequency oji and with a bandwidth T, the bath correlators are 
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Because the qubit is strongly detuned with the cavity, the squeezing terms will only have a negligible effect on the 
qubit. The qubit is thus coupled to a thermal bath with the population ha = sinh r ^ 2 _)_p 2 and the damping rate 7 ^. 
One can then compute the dynamics of (ct^) 
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The qubit excited state relaxes with the Purcell decay rate 73 . 

As expected, the contribution from ha to the decay is therefore largely suppressed for A ^ T and hence pre¬ 
squeezing (i.e. to 7 ^ 0) does not affect the qubit relaxation rate for T ^ A. 


B. Purcell decay in the straddling regime 

To obtain the opposite-signed dispersive shifts required for the TMSS protocol, one of the readout cavities is parked 
in the straddling regime with the qubit, i.e. A < 27r|/ge — /e/|, where fge, fef denote the lowest two qubit transition 
frequencies. The qubit-resonator coupling strength, however, can be chosen to be proportionately weaker: wlO- 

20 MHz instead of the usual 100-200 MHz used in the standard dispersive regime. As a result, the smaller value of 
A is compensated by the smaller value of g, and the Purcell decay rate is no larger than in the usual case (where 
both g and A are larger). At the same time, due to the enhanced contributions from the second excited state, a large 
dispersive shift can still be realized to match the shift induced by the dispersively coupling cavity. Thus, due to the 
value of (//A « 0.1 being similar for both the readout cavities, both will have similar Purcell decay rates. 

For optimal readout conditions, i.e. n = 2%, and assuming the squeezing bandwidth to be sufficiently small (P/A < 
0.1) so that any decay due to extra photons generated by squeezing can be ignored, we estimate the Purcell decay 
for the configuration presented in Figure 4(b) of the main text. For Ej/h = 25 GHz, Ec/h = 0.325 GHz, gij^TT = 
8 MHz, g^l^-K = 15 MHz, we obtain a (single mode) Purcell decay time of 18 /iS (for Cavity 1 at 7.6 GHz) and 11 
(for Cavity 2 at 7.9 GHz), corresponding to (calculated) qubit frequencies /ge = 7.72 GHz, /g/ = 7.35 GHz. Given the 
short measurement times [c.f. Figure 3(a) of main text], there should be no appreciable qubit decay, and concomitant 
loss of readout fidelity, during the measurement. 

We note that it is also possible to implement the straddling regime with strongly anharmonic qubits (such as the 
flux qubits or fluxonium) having higher anharmonicities of few GHz. This allows greater leeway with respect to 
qubit-resonator detunings and dispersive shift profiles. This was exploited in the experiment reported in [S5] where 
the straddling regime was realized with a flux qubit. Despite the very large dispersive shifts (~ 80 MHz), the authors 
of that paper have reported no degradation of the qubit relaxation time due to resonator-induced decay. 


IV. EFFECT OF LOSSES 

A. Photon loss 

We hrst study the effects of photon losses from the cavity output fields by adding a beam splitter at the output 
of each cavity j = 1 , 2 ; these beam splitters combine the cavity output field (which is squeezed) with vacuum 

noise (operator Cin,g). Letting ij denote the transmittance, the output field Aout,j of each beam splitter is Aontj = 
^/V&outj + Vl — rjCinj- The output field ^Vout.-i- leaving the beam splitters then reads 

Jkout,-t- — '\/v^out,+ T 's /1 g yin ,+, (S25) 

where yin,+ = 1 ~ Dn.i + 2 ~ Cin, 2 ) is the joint quadrature of the vacuum fluctuations that has zero mean 

and correlations {yin,+ {t)yin,+ (t')) = S{t — t'). The vacuum correlations of the joint quadrature actually correspond 
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to the output signal correlations without squeezing. The measurement operator after the beam splitters is defined as 
dtyout,+ (i)- From the mean and two-time correlations, the measurement signal and noise read 

= ^/rjM{T), 

-I- (1 — r7)(M^(T))r=0 

The SNR after the beam splitters is finally equal to 
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The general result Eq. (S28) shows that to fully exploit the properties of squeezed states in quantum metrology, 
photon loss has to be made minimal; as mentioned in the main text, this is true for any scheme which attempts to 
exploit squeezing. As the SNR is exponentially enhanced in the absence of photon loss, SNR(r) = e’’SNR(r = 0), we 
get 


SNR,,(r) 


SNR(r) 


(S29) 


For non-zero photon loss, the SNR saturates as a function of r (as opposed to increasing indefinitely). Practically 
speaking, for rj < 1, the maximum useful amount of squeezing corresponds to ^ T^i further increases in r do not 
significantly improve the SNR. 


B. Internal loss 

Photon losses occurring inside the cavities (internal loss) have a qualitatively similar effect to downstream losses. 
We study the effect of an additional decay channel to a reservoir at zero temperature in each cavity. The damping 
rate to the two-mode squeezed reservoir is noted fj,K and the damping rate to vacuum is noted (1 — /r)K, keeping the 
total damping rate k. Following Sec. I A, we first calculate the Langevin equation of the squeezed QMFS, 

^s{t) = I ) Zs{t) -y^Zs,Ut) - y/(W^Zi„(t), (S30) 

V-Xu^ -2V 

where Zi„ are the joint operators of the vacuum baths. The output field is then calculated using the boundary 
condition Zg^out = v^/dcZs -|- Zg^n. The measurement signal turns out to be equal to 

M^,s(t) = m'Ms(t). (S31) 



FIG. SI. Measurement time to reach a fidelity of 99.99% against squeeze strength without loss (red line) and in the presence 
of photon loss, Eq. (S29) (orange line), and internal loss, Eqs. (S31)-(S32) (indigo fine). 
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Again following Sec. I A, the measurement noise is found to be 

|l + (e^'' - l)/r(l - ^)4cos^( 

The spurious decay channel reduces the measurement signal and gives rise to a transient behavior to the measurement 
noise, thereby altering the corresponding SNR^. The effect of the internal loss on the measurement time is plotted in 
Fig. SI and compared to the case of photon loss from the cavity output fields. 
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